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ABSTRACT 

We derive the baryon mass function for a complete sample of low-redshift clusters and argue that it is an ex- 
cellent proxy for the total mass function if the ratio ft, = Mh/M tot in all clusters is close to its universal value, 
Vlb/Q,M- Under this assumption, the baryon mass function can be used to constrain the amplitude and slope 
of the density fluctuations power spectrum on cluster scales. This method does not use observational determi- 
nations of the total mass and thus bypasses major uncertainties in the traditional analyses based on the X-ray 
temperature function. However, it is sensitive to possible systematic variations of the baryon fraction as a func- 
tion of cluster mass. Adapting a weak dependence fb(M) suggested by observations and numerical simulations 
by Bialek et al., we derive <j% = 0.72 ±0.04 and the shape parameter il^h = 0.13 ±0.07, in good agreement with 
a number of independent methods. We discuss the sensitivity of these values to other cosmological parameters 
and to different assumptions about variations in ft,. 

Subject headings: cosmological parameters — galaxies: clusters: general — surveys — X-rays: galaxies 



1. INTRODUCTION 

The number density of galaxy clusters, as a function of 
mass, is an excellent indicator for the normalization and 
slope of the power spectrum of the density fluctuations in the 
present-day Universe. The power spectrum normalization is 
traditionally expressed in terms of a$, the rats amplitude of 
the density fluctuations on the 8/T 1 Mpc scale. In the CDM 
cosmology, the shape of the present-day power spectrum is 
controlled almost completely by the product £Imk (Bond & 
Efstathiou 1984). Therefore, the cluster mass function is a 
sensitive probe of these cosmological parameters. 

Considerable progress in the theoretical models of structure 
formation in the past decade has provided a reliable machin- 
ery to compute the shape of the power spectrum for any given 
set of cosmological parameters (Seljak & Zaldarriaga 1996, 
Eisenstein & Hu 1998), and to predict the mass function of 
collapsed objects (clusters) given the power spectrum (Jenk- 
ins et al. 2001). Most of the problems with using the cluster 
mass function as a cosmological probe are on the observa- 
tional side. 

Theory operates with the virial masses of clusters. How- 
ever, the total cluster mass within the virial radius, as opposed 
to the mass within a small metric radius, is notoriously dif- 
ficult to measure. The main techniques such as weak lens- 
ing, the X-ray method assuming hydrostatic equilibrium of 
the intracluster gas, and velocity dispersion of cluster galax- 
ies, lead to significant statistical and systematic uncertainties 
at large radii, and often produce contradictory results (see e.g., 
Markevitch & Vikhlinin 1997 and Fischer & Tyson 1997 for 
representative uncertainty estimates for the X-ray hydrostatic 
and weak lensing methods, respectively). It is probably fair to 
say that total masses at the virial radius cannot be determined 
to better than 30^-0% accuracy at present. 

Because of the difficulties with the virial mass measure- 
ment, a direct derivation of the mass function for a large, com- 
plete sample of clusters is not feasible at present. Therefore, 
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the cluster mass is usually related to a more easily measured 
quantity. Consider one of the most frequently used proxies for 
the mass function, the cluster temperature function. This re- 
lies on the theoretically expected tight correlation M oc T 3 / 2 . 
The theoretically expected slope of the M-T correlation is 
reproduced both in numerical simulations (Evrard, Metzler & 
Navarro 1996) and observationally, within a given approach 
to the total mass measurement (Horner, Mushotzky, & Scharf 
1999; Nevalainen, Markevitch & Forman 2000; Finoguenov, 
Reiprich & Bohringer 2001). However, the normalizations of 
the M-T relation derived by different authors span a wide 
range (see Evrard et al. 2002 for a review). Ultimately, the 
question of the normalization of the M-T relation boils down 
to the possibility of total mass measurements with small sys- 
tematic uncertainties, which is at the limits of current obser- 
vational techniques. 

The expected cluster mass functions are very steep and 
therefore any change in the mass scale leads to large variations 
of the derived cosmological parameters. In a recent review, 
Rosati et al. (2002) estimated that 30% uncertainties in the 
normalization of the M- T relation leads to a 20% uncertain- 
ties in (78 and 30% uncertainties in Hm derived from cluster 
evolution. Until the M-T relation is reliably normalized, pre- 
cise measurements of the cosmological power spectrum from 
the cluster temperature function are impossible. Therefore, it 
is desirable to explore new proxies for the cluster mass func- 
tion that would not require absolute measurements of the total 
mass. In this Paper, we propose that such a proxy can be based 
on the easily measured cluster baryon masses. 

Indeed, it is expected that the fraction of baryons in the 
virial mass of clusters equals the average value in the Uni- 
verse — f, = M h /M tot = n b /n M (White et al. 1993). There- 
fore, the baryon fraction is the same for clusters of any mass 
at any redshift at least to a first approximation. This is one 
of the best-proved cluster-related facts, that is justified by ba- 
sic theory, confirmed by numerical simulations (e.g., Frenk et 
al. 1998), and supported observationally (Mohr, Mathiesen & 
Evrard 1999, Allen et al. 2002). Given the universality of the 
baryon fraction, there is a trivial relation between the cumu- 
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lative total and baryon mass functions: 

F h (M b ) = F tot (n M /n b M h ). (i) 

The average baryon density in this equation is fixed by ei- 
ther the Big Bang Nucleosynthesis theory or by the cosmic 
microwave background fluctuations, Q b h 2 = 0.0224 ± 0.0009 
(Spergel et al. 2003). The baryon mass is a convenient proxy 
for the total mass because 1) it is easily derived from the X-ray 
imaging observations (§|3Jl, and 2) the scaling between Mb and 
M tot is expressed in terms of the model parameters and 
therefore does not rely on the absolute total mass measure- 
ments. In this Paper, we apply this method to ROSAT obser- 
vations of a complete flux-limited sample of 52 low-redshift 
clusters selected from the All-Sky Survey data. 

The Paper is organized as follows. In §|2] we describe selec- 
tion of our cluster sample. The measurements of the baryon 
masses are presented in §[3] and the baryon mass function is 
derived in §|4] The theory involved in modeling the baryon 
mass function is reviewed in §[5] In §[6] we derive constraints 
on erg and on the shape parameter, T = fli^h, and discuss the 
dependence of these results on our basic assumption of the 
universality of the baryon fraction. 

Numerical values of the cluster parameters are quoted for 
the Hq = 71 km s" 1 Mpc" 1 , Q,m = 0.3 and A = 0.7 cosmology. 

2. CLUSTER SAMPLE 

Our cluster sample was selected from the HIFLUGCS cata- 
log (Reiprich & Bohringer 2002) of clusters selected from the 
ROSAT All-Sky Survey (RASS) data. This catalog includes 
63 clusters above a flux limit of 2 x 10" 11 erg s" 1 cm" 2 in the 
0.1-2.4 keV band, detected over 2/3 of the sky. 

We have applied further objective selections to facilitate a 
more reliable measurement of the baryon mass function. First, 
we increased the flux limit to 1.4 x 10" 11 erg s" 1 cm" 2 in the 
0.5-2 keV band (or 2. 3 x lO^'ergs" 1 cm" 2 in the 0.1-2.4 keV 
band). Second, we considered only clusters at z > 0.01 be- 
cause at lower redshifts, the virial radius of most objects does 
not fit inside the ROSAT field of view. In total, 52 clusters 
from the HIFLUGCS catalog satisfied our selection criteria 
(Table[T). 

Most of the selected objects were observed by ROSAT 
in one or more pointed observations, so high-quality X-ray 
data are available. Several more clusters are so bright that 
their baryon mass can be measured from the shallower RASS 
data. For only 5 clusters, A576, A1736, A2065, A3581, and 
EX0422, the existing X-ray data do not permit a direct baryon 
mass measurement at the virial radius, so we used the Mb —L x 
correlation to estimate their masses. 

3. BARYON MASS MEASUREMENT 
3.1. ROSAT Data Reduction and Gas Mass Measurements 

The data were prepared in a standard manner. For pointed 
observations, we used S. Snowden's software (Snowden et al. 
1994) to produce flat-fielded images with all non-sky back- 
ground components removed. For the RASS data, we used the 
archive-supplied exposure maps. Experiments with different 
energy bands show that the 0.7-2 keV band has the optimal 
ratio of cluster to background surface brightness. Our further 
analysis is performed in this band. 

The measurement of the gas mass profile from the X-ray 
imaging data was performed almost identically to the ap- 
proach described in Vikhlinin, Forman & Jones (1999). We 
first extracted the azimuthally-averaged surface brightness 



profile centered on the apparent centroid of the cluster. The 
background level was estimated by fitting the profile at large 
radii, 0.3r v < r < 1.5r v where r v is the virial radius estimated 
from the X-ray temperature 3 , with the power law plus con- 
stant model. The power law here represents the cluster bright- 
ness and the constant corresponds to the background. The 
background intensity was subtracted from the profile assum- 
ing that it is constant across the image, 

Assuming that the clusters are spherically symmetric, the 
gas density profiles were derived from the surface brightness 
profiles using direct deprojection (Fabian et. al. 1981). The 
deprojection method yields the radial profile of the X-ray vol- 
ume emissivity which is easily converted to the gas density as 
n e = Ke 1 / 1 , p g = \ .\5n e m p , where the numerical value for K 
can be found for any plasma temperature using the MEKAL 
code (Mewe et al. 1985), and the numerical value in the sec- 
ond relation corresponds to the standard mix of fully ionized 
hydrogen and helium. 

For three clusters, A 1060, A 1644, and NGC1550, the 
ROSAT surface brightness profile were too noisy for direct 
deprojection. We resorted to the usual /3-model analysis in 
these three cases. 

The intracluster gas which we directly observe in X-rays is 
the dominant baryonic component of clusters. We will use the 
cluster masses corresponding to the density contrast i5 = 324 
relative to the mean density of the Universe at the redshift 
of observation (see §[3]below). Therefore, we derive the gas 
masses from the following equation: 3M g (< r)/(4irr 3 (pb)) = 
324, where M g ( <r) is the gas mass profile as derived from the 
X-ray data. We adopt the average baryon density as given by 
the WMAP observations of the cosmic microwave background 
fluctuations, (p b ) = 6.23 (1 + z) 3 M Q kpc" 3 (or Q b h 2 = 0.0224, 
Spergel et al. 2003) which is also close to the Big Bang Nu- 
cleosynthesis value of ilb h 2 = 0.020 (Buries, Nollett & Turner 
2001). The observed X-ray surface brightness profiles permit 
a direct measurement of the gas mass at the radius correspond- 
ing to S = 324. 

The measured gas masses M g 324 are listed in Tabled It 
is useful to provide simple relations which allow one to scale 
our results to other values of the mean overdensity, Hubble 
constant etc. Empirically, we find that on average, the radial 
dependence of mass is M(9) oc 9, which corresponds to the 
following scaling as a function of overdensity 

MsocfT ' 5 . (2) 

Within a fixed aperture, the X-ray derived baryon mass scales 
as /i~ 5//2 . There is an additional scaling because Ms is mea- 
sured at a fixed overdensity 5 rather than within a fixed an- 
gular radius. The angular radius corresponding to the over- 
density 5 is found from the equation M(6s)/(0lh~ 3 ) = const. 
Given the empirical average mass profile, M(9) oc 9, and the 
/z-scaling of M(6), 9s should scale as h° 25 , and therefore 

Maoc/r 2 - 25 . (3) 

The only effect of changing the average baryon density in the 
Universe is to change the apparent mean overdensity, and so 
one can use eq. (|2jl to convert our masses to 8 = 324 for a 
different value of il b h 2 : 

Ms^iflbh 2 )- 05 . (4) 

3 We used the relation provided by Evrard, Metzler & Navarro (1996), 
r, = 1.95ft" 1 MpcCr/lOkeV) 1 / 2 
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Table 1 
Cluster sample 



Object 


z 




fx" 


L x 






r «,324 


£-opt C 






(keV) 




10 44 erg s- 1 


10 14 M Q 


1O 14 M 


(Mpc) 


1O 13 L 


2A0335 


0.0349 


3.6 


6.8 


1.8 


0.59 ±0.07 


0.78 ±0.10 


1.85 




A85 


0.0556 


6.9 


4.3 


3.0 


1.11 ±0.13 


1.35 ±0.16 


2.24 


0.41 ±0.06 


A119 


0.0442 


5.6 


2.4 


1.0 


1.00 ±0.10 


1.27±0.13 


2.19 




A262 


0.0155 


2.2 


4.6 


0.2 


0.20 ±0.04 


0.31 ±0.06 


1.31 


0.09 ±0.03 


A399 


0.0715 


6.8 


1.7 


2.0 


1.21 ±0.31 


1.62 ±0.38 


2.27 




A400 


0.0238 


2.3 


1.6 


0.2 


0.25 ±0.02 


0.38 ±0.04 


1.41 


0.11 ±0.03 


A401 


0.0739 


8.0 


3.1 


3.9 


2.34 ±0.26 


2.82 ±0.31 


2.82 


0.47 ±0.18 


A478 


0.0882 


6.9 


4.5 


8.2 


2.16±0.24 


2.59 ±0.29 


2.71 


0.36 ±0.09 


A496 


0.0328 


4.7 


5.4 


1.3 


0.61 ±0.08 


0.80 ±0.11 


1.87 




A576 


0.0389 


4.0 


1.6 


0.5 




0.46 ±0.14* 






A754 


0.0542 


9.8 


4.3 


2.9 


1.58±0.16 


1.90 ±0.20 


2.52 




A1060 


0.0137 


3.2 


6.1 


0.2 


0.20 ±0.04 


0.29 ±0.06 


1.31 




A 1367 


0.0214 


3.7 


4.7 


0.5 


0.41 ±0.07 


0.54 ±0.09 


1.66 


0.14 ±0.04 


A1644 


0.0474 


4.7 


2.2 


1.1 


0.92 ±0.21 


1.19±0.25 


2.11 




A1651 


0.0860 


6.2 


1.5 


2.6 


1.14±0.17 


1.42 ±0.21 


2.19 




A1656 


0.0231 


8.1 


18.2 


2.1 


1.54±0.18 


1.91 ±0.23 


2.57 


0.51 ±0.11 


A1736 


0.0461 


3.6 


1.7 


0.8 




0.66 ±0.20* 






A 1795 


0.0622 


6.0 


3.9 


3.4 


0.93 ±0.08 


1.15 ±0.11 


2.10 


0.25 ± 0.04 


A2029 


0.0766 


9.1 


4.3 


5.7 


1.86 ±0.24 


2.30 ±0.29 


2.60 


0.58±0.12 


A2052 


0.0353 


3.4 


3.0 


0.8 


0.31 ±0.05 


0.42 ±0.07 


1.49 


0.13 ±0.03 


A2063 


0.0355 


3.7 


2.5 


0.7 


0.44 ±0.07 


0.60 ±0.09 


1.68 


0.17 ±0.03 


A2065 


0.0726 


5.4 


1.4 


1.7 




1.13±0.34* 






A2142 


0.0894 


9.3 


3.9 


7.3 


3.13 ±0.23 


3.73 ±0.28 


3.06 


0.56 ±0.09 


A2147 


0.0353 


4.4 


3.1 


0.9 


1.01 ±0.18 


1.33 ±0.22 


2.21 




A2163 


0.2010 


11.5 


1.4 


14.6 


5.58 ±0.74 


6.41 ±0.85 


3.37 




A2199 


0.0299 


4.8 


6.5 


1.3 


0.57 ±0.04 


0.75 ±0.06 


1.84 


0.24 ±0.04 


A2204 


0.1523 


7.1 


1.6 


9.4 


2.07 ±0.43 


2.51 ±0.49 


2.52 




A2256 


0.0581 


7.3 


3.7 


2.8 


1.48±0.12 


1.80±0.15 


2.45 


0.54 ±0.06 


A2589 


0.0416 


3.7 


1.6 


0.6 


0.29 ±0.04 


0.42 ± 0.06 


1.46 




A2634 


0.0309 


3.3 


1.5 


0.3 


0.31 ±0.04 


0.42 ±0.05 


1.50 


0.19 ±0.03 


A2657 


0.0400 


3.7 


1.7 


0.6 


0.40 ±0.05 


0.54 ±0.06 


1.62 




A3112 


0.0750 


5.3 


1.9 


2.5 


0.81 ±0.13 


1.03±0.16 


1.98 




A3158 


0.0591 


5.8 


2.3 


1.8 


0.90±0.15 


1.17±0.19 


2.09 




A3266 


0.0589 


8.0 


3.5 


2.8 


1.71 ±0.17 


2.05 ±0.20 


2.57 




A3376 


0.0455 


4.0 


1.4 


0.6 


0.41 ±0.06 


0.54 ±0.07 


1.63 




A3391 


0.0514 


5.7 


1.5 


0.9 


0.93 ±0.14 


1.24±0.18 


2.12 




A3395 


0.0506 


5.0 


1.9 


1.1 


1.24 ±0.20 


1.59 ±0.25 


2.33 




A3526 


0.0112 


3.4 


15.7 


0.4 


0.31 ±0.05 


0.44 ±0.07 


1.51 




A3558 


0.0480 


5.5 


4.0 


2.1 


1.34±0.12 


1.66±0.16 


2.39 




A3562 


0.0490 


5.2 


1.8 


1.0 


0.92 ±0.09 


1.17±0.13 


2.12 




A3571 


0.0391 


7.2 


7.5 


2.5 


1.21 ±0.18 


1.50 ±0.22 


2.33 




A3581 


0.0214 


1.8 


2.0 


0.2 




0.23 ±0.07* 






A3667 


0.0556 


7.0 


4.6 


3.2 


2.50 ±0.24 


2.99 ±0.29 


2.93 




A4038 


0.0292 


3.3 


3.6 


0.7 


0.31 ±0.05 


0.44 ±0.07 


1.50 




A4059 


0.0460 


4.4 


2.0 


0.9 


0.45 ±0.07 


0.62 ±0.09 


1.67 




EXO0422 0.0390 


2.9 


1.9 


0.6 




0.53±0.16* 






Hydra A 


0.0522 


3.6 


2.9 


1.8 


0.60 ±0.07 


0.77 ±0.09 


1.83 




MKW3s 


0.0453 


3.7 


2.1 


1.0 


0.40 ±0.06 


0.53 ±0.08 


1.60 




MKW4 


0.0210 


1.7 


1.4 


0.1 


0.14 ±0.01 


0.21 ±0.03 


1.15 


0.05 ±0.02 


NGC1550 


0.0120 


1.1 


2.5 


0.1 


0.05 ±0.01 


0.08 ±0.02 


0.84 




NGC507 


0.0155 


1.3 


1.5 


0.1 


0.08 ±0.01 


0.13 ±0.02 


0.97 




S 1101 


0.0580 


2.6 


1.5 


1.1 


0.29 ±0.04 


0.39 ±0.06 


1.43 





" Cluster temperatures adapted from Markevitch et al. (1998, 1999), White et al. (2000), Fukazawa et al. (1998), David et al. 
(1993), Arnaud et al. (2001), Blanton et al. (2001), Horner et al. (1999), Allen & Fabian (1998), Finoguenov et al. (2001), 
Kaastra et al. (2001), Sun et al. (2003) 

* X-ray flux in 0.5-2.0 keV, in 10~"erg s -1 cnT 2 units. 

c Optical luminosity inside r g ^ 2 4 in the V band, adapted from Arnaud et al. (1992), Hradecky et al. (2000), Beers et al. (1984). 

* — baryon mass determined from M4324-LA correlation. 
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FIG. 1. — Correlation of the optical luminosity of the cluster in the V 
band and the gas mass (M g .324). The original optical measurements were 
extrapolated to r g 324 as explained in the text. The solid line shows the best- 
fit power law (eq.jSj 



3.2. Stellar Mass in Clusters 

The stellar material in galaxies contributes a small but non- 
negligible fraction of the cluster baryon mass. Not all clusters 
in our sample have high-quality optical observations. There- 
fore, we first use the published data for a subsample of our 
clusters to establish a correlation of the optical luminosity 
with the gas mass and then use this correlation to estimate 
the stellar contribution in all clusters. 

The optical luminosities for 16 clusters were compiled from 
the papers by Arnaud et al. (1992), Hradecky et al. (2000), and 
Beers et al. (1984). The optical luminosities in these papers 
are quoted at different metric radii in the V band, and so we 
converted them to the gas overdensity radius ryi\ assuming 
that the light follows the King profile with r c = 180 kpc (Bah- 
call 1975). The resulting luminosities within 7-324 show a very 
tight correlation with the gas mass (Fig.[0 which is fitted by 
the following relation 



L opt = 319L 



M„ 



324 



0.71 



(5) 



We further assume the same mass-to-light ratio for the stars, 
M*/L opt = 6Mq/Lq in the V band, which follows from the 
stellar evolution models for the average population of ellipti- 
cal and spiral galaxies in clusters (Arnaud et. al. 1992). Under 
these assumptions, equation can be used to estimate the to- 
tal baryon (i.e. gas plus galactic stars) mass of a cluster from 
the X-ray measured gas mass: 



M fc ,324 



M, 



= 1.100 + 0.045 



S,324 



M, 



S.324 



10 15 M^ 



-0.5 



(6) 



Therefore, galactic stars contribute about 15% of the total 
baryon mass in massive clusters, and this fraction slightly 
increases in the less massive systems. This estimate of the 
stellar contribution is consistent with the considerations pre- 
sented in Fukugita, Hogan & Peebles (1998). The trend in 
the stellar-to-gas ratio is similar to that derived by Lin, Mohr 
& Stanford (2003) using the Zf-band photometry from the 



A, 



Kf 



10 1 - 5 



I0 4 



to 44 

L/(l+:) 2 - 26 ,erg S - 



FlG. 2. — Cluster baryon mass as a function of the total X-ray luminosity in 
the 0.5-2 keV band. The solid line shows the best fit power law and dashed 
lines indicate the 1<t scatter around the mean relation. 



2MASS survey, although the overall normalization of M*/M g 
is smaller in that paper. 

We ignore any contribution from intergalactic stars to the 
cluster baryon mass. At present this component cannot be 
reliably estimated. However, all available studies (e.g., Feld- 
meier et al. 2002) suggest that the intergalactic stars contribute 
only a small fraction to the total stellar mass of the cluster, and 
therefore only a few percent at most to the total baryon mass. 

3.3. Baryon Mass Error Budget 

The main contribution to the mass measurement uncertainty 
is the Poisson noise in the X-ray surface brightness profiles. 
This noise component is straightforward to compute in each 
case. Typically, the statistical uncertainty in the baryon mass 
is about 20%. All other uncertainties are systematic and we 
can only estimate their contribution. 

We measured the baryon masses under the assumption that 
the clusters are spherical. Any deviations from the spherical 
symmetry result in a scatter in our mass determination. We 
estimate this scatter to be 7% based on the azimuthal scat- 
ter of the gas mass profiles (Vikhlinin et al. 1999) and from 
the results of numerical simulations by Mathiesen, Evrard & 
Mohr (1999). 

The stellar mass was estimated using a number of assump- 
tions. First, we assumed a standard mass-to-light ratio for 
the galactic stars. Second, the total light was estimated from 
the extrapolation of the King profile to large radii. Third, for 
most clusters, the optical luminosity was estimated from the 
L opt -Mg relation which has some scatter (Fig.[Q. We con- 
servatively estimate that the stellar mass is measured with a 
30% uncertainty. The stellar mass is only 15-20% of the gas 
mass, therefore the contribution to the uncertainty in the total 
baryon mass is small (5-7%). 

The final uncertainties on the baryon mass were obtained 
by adding the above three components in quadrature. The 
dominant contribution is still the statistical noise in the X-ray 
data. 

3.4. Correlation of the X-ray Luminosity and the Baryon 

Mass 
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FIG. 3.— Deviations of the measured luminosity from the mean U b -L x FIG. 4.— Sample volume as a function of the cluster baryon mass, 

relation. Solid line shows the Gaussian fit. 



The key relation for estimating the mass functions is the 
volume of the survey as a function of mass. For a flux-limited 
sample like ours, the volume is determined by the scaling re- 
lation between the mass and the X-ray luminosity. Such a 
relation is also useful for estimating the masses in those cases 
where direct mass measurements are impossible. 

Our sample was selected using the total (i.e. without ex- 
cluding either substructures or cooling flow regions) X-ray 
fluxes in the 0.5-2 keV band. Therefore, the Mt,-L x correla- 
tion should be with the total luminosity in this energy band. 
The corresponding measurements are listed in Table and 
plotted in Fig.|2] The average correlation is well 

fitted by the power law 4 



M/,,324 = 8.5 x 10 1J M Q x (L t /10 44 erg s" 1 ) 



K0.78±0.04 



(7) 



where the best-fit parameters were found using the bisector 
method of Akritas & Bershady (1996) and the uncertainty in 
the slope was estimated using bootstrap resampling. 

There is a significant scatter in the - L x correlation 

which can be partly explained by the contribution of cooling 
flow regions and substructures to the total cluster flux. As Fig- 
ure shows, the scatter around the mean relation is well ap- 
proximated by the log-normal distribution with an rms scatter 
of cigM = 0.19 in mass. Part of this scatter is due to the mass 
measurement errors. Subtracting this contribution in quadra- 
ture, we find that the intrinsic log scatter in the Mb ,324 _ Lx 
correlation is o\ % m = 0.17 in mass for a given luminosity or 
cigz. = 0.15 in luminosity for a given mass. 

4. OBSERVED BARYON MASS FUNCTION 

The volume of a flux-limited survey for the cluster of lumi- 
nosity L is obtained by integration of the cosmological depen- 
dence of comoving volume, dV / dz, over the redshift interval 
Zmin-Zmax, where z m i n is the lower redshift cut (0.01 in our 
case), and z max is the maximum redshift for the cluster to ex- 

4 The M b -L x relation evolves strongly to 7. ~ 0.5 (Vikhlinin et al. 2002). 
To correct for this effect, we divided the luminosities by the observed rate of 
the evolution (1 +z) 2 ' 26 even though this correction is small in our redshift 
range. 



ceed the flux limit, f\\ m : 
/ lim (0.5-2): 



L(0.5(l+z max )-2(l+z max )) 
47r^(z max ) 



(8) 



where <iz.(z max ) is the luminosity distance as a function of red- 
shift. The solution of this equation with respect to z max leads 
to the survey volume as a function of cluster luminosity, V(L). 

Given the probability distribution p(L\M) of the cluster lu- 
minosity for the given mass, the survey volume as a function 
of mass is computed as V(M) = J p(L\M)V(L)dL. Taking into 
account the average relation between M& and L (eq.0 and the 
log-normal scatter around this relation, we obtain: 



V(M): 



1 



(27T) 1 / 2 (TlgM 

(lgL/W-lgD 21 



exp 



2cr lgM 



(9) 



V(L)dlgL, 



where L/(M) is given by the inverse of eq. Q. The sample 
volume as a function of the cluster baryon mass is shown in 
Fig.E 

We will use the cumulative representation of the mass func- 
tion which is estimated from observations as 

f(M)=N(>M)= 777777 j (10) 



M,>M 



V(Md ■ 



where V(M,) is given by eq. |9). The resulting baryon mass 
function is shown in Fig. [5] 

There are two contributions to the uncertainty in the f{M) 
measurement — the Poisson noise which is trivially estimated 
as 



(Apf(M)) 2 = Y, 



1 



M,>M 



V 2 (M,)' 



(11) 



and the component related to the mass measurement un- 
certainties. To take both components into account, we 
performed the Maximum-Likelihood fits of the measured 
f(M) with a two-parameter Schechter function df/dM = 
A(M/M ) a exp(-M/M ) (Schechter 1976). The likelihood 
function was constructed to include the scatter due to the mass 
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measurement errors (for details, see Voevodkin, Vikhlinin & 
Pavlinsky 2002 and Vikhlinin et al. 2003). The range of sta- 
tistically acceptable parameters A, a, and Mo defines a band 
of the cumulative mass functions consistent with the data at 
the given confidence level. The resulting 68% uncertainties 
in the f(M) measurement are shown as a grey band in Fig. [5] 
For reasons explained below (§|6jl, we use only clusters with 
baryon mass Mb > 4 x 10 13 M Q for the theoretical modeling. 
Therefore, the Schechter function fits were performed only in 
this mass range. 

In addition to the statistical uncertainty in the mass function 
there is an additional component because of the cosmic vari- 
ance in a survey of a finite volume. Hu & Kravtsov (2003) 
calibrate the ratio of the cosmic variance and Poisson noise 
as a function of cluster mass. Using their results, we find 
that the cosmic variance is negligible in our mass range for 
a% ~ 0.7-0.8 indicated by our data (see below). 

5. REVIEW OF THEORY 

Recent numerical simulations of structure formation in the 
Universe provide accurate analytic expressions for the cluster 
total mass function (Jenkins et al. 2001). Jenkins et al. have 
demonstrated that the mass functions in the CDM cosmology 
has a universal form when expressed as a function of the linear 
rms density fluctuations on the mass scale M, a(M): 



M dN(>M) 
Jp)d\oga(M) 



■■ A exp (- 1 log a' 1 + b | c ) (12) 



All dependencies of the mass function on cosmological pa- 
rameters enter through a(M). The function a(M) is equivalent 
to the power spectrum of linear density perturbations (Pee- 
bles 1981). We will assume that the power spectrum, P(k), is 
the product of the primordial inflationary spectrum k" and the 
CDM transfer function: P(k) ~ k"T 2 (k). The transfer function 
depends on the parameters Hm, h, and f^. We used analytic 
approximations to T(k) provided by Eisenstein & Hu (1998). 
Therefore, the theoretical model for the mass function is de- 
fined by parameters CIm, f2&, h, n, and also the normalization 
of the power spectrum, a%. 



The numerical values of coefficients A, b, and c in dl2l de- 
pend on the definition of the cluster mass. The most con- 
venient definition, from the observational point of view, is 
the mass corresponding to a given spherical overdensity, 6, 
relative to the mean density at the cluster redshift. Jenkins 
et al. considered two values of 8, 180 and 324 (the latter 
value corresponds to the virial radii in the ACDM model with 
Qm = 0.3). We measured baryon masses for S = 324. For 
this overdensity, Jenkins et al. obtain A = 0.316, b = 0.67, and 
c = 3.82. 

If the baryon fraction in clusters is indeed universal, the 
baryon and total overdensities must be equal, 5 = 5b, which 
allows one to use the measured baryon mass Mb32A as a proxy 
for the total mass, M324. Recall that the model for the baryon 
and total mass function are related via equation Q. 

Using the machinery described above, we can compute the 
model for the baryon mass function for any set of cosmo- 
logical parameters 5 and then fit it to the observations, self- 
consistently rescaling the observed mass function to the cos- 
mological parameters being tried. In the rest of this section, 
we consider separately the effect of each cosmological param- 
eter on the Fb(Mb) model and on scaling of the observed mass 
function. We also consider the effect of possible deviations of 
Mb /M tot from the universal value. 

5.1. Parameter Qm 

Parameter Q,m, mainly through the product ^ImK deter- 
mines the shape of the transfer function T(k) and therefore 
the slope of the mass function. This is a strong effect which 
allows us to derive VIm from the cluster mass function. 

In addition, Q,m affects the growth of density perturbations. 
This effect should be taken into account, because the observa- 
tional mass function is derived in a redshift interval of finite 
width. Therefore, we computed the mass function models at 
several redshifts and weighted them with the survey volume 
for the given mass. 

The effect of VIm on the observed mass function via the 
mass and volume calculations is very weak and can be ig- 
nored. The growth factor also depends on A, but very weakly 
at low redshift, and therefore we fixed A = 0.7. 




FIG. 5. — Local baryon mass function. The gray band shows a 68% uncer- 
tainty interval derived from the Maximum Likelihood analysis including the 
mass measurement errors. The solid line corresponds to the prediction of our 
best-fit cosmological model (Fig.lol. 



5.2. Power Spectrum Normalization a% 

The normalization of the mass function model is exponen- 
tially sensitive to erg, which allows a precise determination of 
this parameter from the observed cluster abundance. 

The baryon mass measurements do not depend on erg, so 
there is no associated scaling of the observed mass function. 

5.3. Hubble Constant h 

The derived baryon masses are very sensitive to the value of 
h, Mb.324 oc h~ 2 25 (eq. and so the observed mass function 
must be scaled accordingly. In the mass function model, the 
Hubble constant enters through the product fi^/i (see above) 
and also changes the mass scale, M oc h~ l (see e.g. Jenkins et 
al.). The normalization of the model mass function scales as 
/1 3 , and, obviously, the observed mass function scales identi- 
cally. 

Although the Hubble constant measurements now converge 
to h = 0.71 (Freedman et al. 2001, Spergel et al. 2003), our 
method is quite sensitive to h. To determine the appropriate 

5 We used the code kindly provided by A. Jenkins. The code was slightly 
modified to include the power spectrum model by Eisenstein & Hu (1998). 
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scalings, we varied h in the range from 0.5 to 0.8. We use 
h = 0.71 for a baseline model. 

5.4. Average Baryon Density VLbh 2 

The average baryon density in the Universe is known quite 
accurately, Q, b h 2 = 0.0224 ± 0.0009 (Spergel et al. 2003, see 
also Buries et al. 2001). This para meter enters the measure- 
ment of the baryon masses Mb .324 (§ 13.11 eq. ®), the scaling 
between the models for the total and baryon mass functions 
(eq. 0, and, weakly, the shape of the transfer function. Vari- 
ations of Vlbh 2 within the above uncertainty range result in 
negligible changes in either the observed or theoretical mass 
function. We, however, allowed f2/,/z 2 to vary in a broader 
range, from 0.018 to 0.025. If not stated explicitly, all the 
results below are reported for Vtbh 2 = 0.0224. 

5.5. Primordial Power Spectrum Slope n 

The observed mass function is independent of the power 
spectrum. For the model mass function, the effect of n 
is indistinguishable from that of the product QmIi because 
our entire mass range corresponds to a narrow (factor of 2) 
wavenumber range, where the power spectrum of density per- 
turbations is indistinguishable from a power law. Our param- 
eter constraints below were generally obtained for n = 1 , ex- 
pected in the inflationary models (Starobinskii 1982, Guth & 
Pi 1982, Hawking 1982) and favored by the recent CMB ob- 
servations (Spergel et al. 2003). The n-dependence of the de- 
rived parameter values was studied by varying n in the range 
from 0.85 to 1.15. 

5.6. Non-Universality of the Baryon Fraction T 

Our method depends on the assumption that the baryon 
fraction in clusters is universal, Mb / M tot = flb/ ^m- Any large 
deviations from the universality may have serious observa- 
tional and methodological implications, and therefore require 
special attention. 

Consider first how the baryon mass measurements have to 
be scaled, if the baryon fraction is non-universal. Let us as- 
sume that at large radii, baryons still follow the dark matter, 
i.e. Mi,/M tot = T ttb/ilM where T is a constant for each clus- 
ter. The first obvious correction is to multiply the baryon mass 
by T . However, if T 1, the baryon and total overdensities 
at each radius are different, 5b = T5. We need the baryon mass 
which corresponds to the total overdensity 5 = 324, and there- 
fore the baryon mass measurements must be corrected by an 
additional factor of T"° 5 (cf. eq.[2j. Overall, the baryon mass 
measurement must be scaled by a factor of T . 

Consider now the effect of possible variants for the non- 
universality of the cluster baryon fraction. Most of the devi- 
ations of T from unity can be factorized into the following 
three possibilities: 

1) T y 1, but is the same for all clusters. Such a devi- 
ation accounts for any systematic under- or overabundance 
of baryons in clusters. Such a deviation is often observed in 
numerical experiments (Mathiesen et al. 1999 and references 
therein). Any systematic bias in the X-ray based gas mass 
measurement is also indistinguishable from T ^ 1 . For ex- 
ample, Mathiesen et al. show that a systematic 5-10% over- 
estimation of the gas mass is possible because of deviations 
from spherical symmetry. Overall, 15% variations of T seem 
to include the range for all known effects. We studied the de- 
pendence of our results on such a systematic change of T by 
repeating all the analyses for T = 0.85, 1.0, and 1.15. 



2) T = 1 on average, but there is some cluster-to-cluster 
scatter. This can be represented by a convolution of the baryon 
mass function model with the appropriate kernel. We find, 
that the effect on the derived parameters is negligible as long 
as the scatter of T remains small (< 15%), which is supported 
observationally (Mohret al. 1999). 

3) T is a function of the cluster mass and T — > 1 in mas- 
sive clusters. This type of variation of the baryon fraction is 
expected, e.g., in the models which involve a significant pre- 
heating of the intergalactic medium. As a toy model for T(M), 
we considered the parameterization 

T(M fe )=l-^,M* = O.7xlO 13 M , (13) 
Mb 

which describes the results of numerical simulations by 
Bialek et al. (2001) where the level of preheating was adjusted 
so that the simulated clusters reproduce the observed L x — T re- 
lation. This relation also is reasonably consistent with recent 
direct measurements of the baryon fraction in a small number 
of clusters spanning a temperature range of 1 .5-10 keV (Pratt 
& Arnaud 2002, 2003; Sun et al. 2003, Allen et al. 2002). 

Equation 1131 was used as our baseline model in deriving 
the parameter constraints, but we also considered all of the 
above possibilities for deviations of T from unity. 

The results of this section can be briefly summarized as fol- 
lows. Given a set of parameters VIm, h, 0%, «, fib, we compute 
the linear power spectrum of density perturbations (Eisen- 
stein & Hu 1998), then convert it to the model for the total 
mass function (Jenkins et al. 2001) and apply to the observed 
baryon mass function scaled along the mass axis as follows 
(here, ujb = Slbh 2 ) 

M/r i = T^Mb^ ( —V 2 ' 25 ( (14) 
Ub \0.1lJ V - 0224 / 

6. RESULTS 

6.1. Application of the Models to the Data 

In the previous section, we discussed how to compute a 
model for the baryon mass function for a given set of CDM 
model parameters. Parameter constraints can be obtained by 
fitting such models to the observed mass function using the 
Maximum Likelihood (ML) method. Note, however, that an 
exact computation of the likelihood in our case is rather cum- 
bersome because of the presence of the mass measurement 
uncertainties in addition to the Poisson noise (§0}. To avoid 
the computational overhead, we used a simpler technique to 
derive the cosmological parameter constraints. We use 68% 
and 95% confidence bands for the mass function obtained in 
§0]using an ML fit to the Schechter function. A mass function 
model which is entirely within a 68% confidence band in the 
interesting mass range is deemed acceptable at the 68% confi- 
dence level (and likewise for the 95% confidence level). Fur- 
ther comparison of this and ML methods is given in Vikhlinin 
et al. (2003). 

This method results in slightly more conservative confi- 
dence intervals than a direct ML fitting. An advantage of our 
approach over ML fits is its direct relation to the goodness- 
of-fit information. By design, all acceptable mass function 
models actually fit the data. On the contrary, the goodness-of- 
fit information is generally unavailable with the ML approach 
and sometimes it is possible to obtain artificially narrow confi- 
dence intervals on parameters only because the best fit model 
provides a poor fit to the data. 
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FIG. 6. — Confidence intervals (68% and 9 5%) for the baseline model 
(h = 0.71, n= 1, tt b h 2 = 0.0224, and T from eq. fTIl ). 



For parameter constraints, we considered the observed 
baryon mass function in a restricted mass range, 4 x 
10 13 M Q <M b <3x 10 I4 M Q . The choice of the lower limit 
is motivated by the possibility that the baryon fraction drops 
systematically in the less massive clusters (see above). Also, 
the number of clusters in our sample with Mi, < 4 x 10 13 Mq is 
small and so the effects of incompleteness may be important. 
The upper limit of the mass interval was chosen to exclude 
the two most massive and potentially peculiar clusters A2142 
andA2163(Fig.|5}. 

As was discussed in §|5] the most sensitive constraints from 
the cluster baryon mass function are for parameters a% and 
VLmIi. Other parameters, h, n, flbh 2 , cannot be constrained 
from such an analysis. Their variation causes simple scalings 
of the constraints on erg and Q,Mh. Below, we will describe 
constraints on <jg and Vl M h for the baseline model, h = 0.71, 
n = 1, VLj,h 2 = 0.0224, and provide scalings of these constraints 
corresponding to variations of the latter three parameters as 
well as deviations of T(M) from the adopted model. We con- 
sider separately the <r%-Q,Mh constraints from fitting the shape 
of the mass function in the 4 x 1O 13 M <Mj<3x 10 I4 M© 
mass range and from using only the normalization of the mass 
function at M b = 1 .0 x 10 14 M Q . 

6.2. Fitting the Shape of the Baryon Mass Function 

A simultaneous determination of the parameters <7g and 
Vl M h is possible from fitting the mass function in a broad Ar- 
range because the slope of the mass function is sensitive to 
fliwh and its normalization — to <7g. In Fig.|6] we show the 
68% and 95% confidence regions for our baseline model. The 
best fit is obtained for VLmIi = 0.13 and erg = 0.72 (the corre- 
sponding mass function model is shown by the solid line in 
Fig. [5}. The 68% uncertainty interval on fti^h includes the 
'Cosmic Concordance' value of VLm = 0.3, but the confidence 
regions are skewed towards lower values of This result 
is similar to a conclusion from the cluster temperature func- 
tion studies that the cluster mass function is shallower than 
predicted in CDM (Henry 2000). 

Parameter h: Consider now how the erg - fti^h constraints 
depend on the exact value of h. The effect of the Hubble con- 
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FIG. 7. — The effect of variations in the baryon fraction on the cos mo- 
logical constraints. The solid contour corresponds to T given by eq. (H) 
(baseline model), the dotted contour is obtained for t he case of a 10% cluster- 
to-cluster scatter in T around the mean relation 1131 . and the dashed contour 
corresponds to T = 1 . 



stant on the shape of the power spectrum is important only 
through the product il^- Therefore, the effect of h is to shift 
the confidence levels in Fig.[6]along the erg axis. Experiment- 
ing with different values of h, we find that this shift is 



Acr 8 = 0.31A/i. 



(15) 



Note that the eg constraints derived from the cluster tem- 
perature function are independent of h because both the mass 
corresponding to a given T and the independent variable of 
the model mass function scale as h~ l . In our case, a different 
scaling of the baryon mass measurement with h (Mb oc /T 2,25 , 
see §|5|l leads to a dependence of the derived erg on h. 

Parameter n: The shape of the cluster mass function is de- 
fined by the slope of the linear power spectrum of density per- 
turbations in the narrow wavenumber range around k ~ 0.1/z, 
where P(k) is accurately approximated by a power law. This 
means that we cannot disentangle the slope of the primordial 
(inflationary) power spectrum n from the slope of the transfer 
function on cluster scales. For example, for a higher value 
of n we need a flatter transfer function, i.e., a lower value of 
VlMh. Therefore, our erg - Vti^h constraints computed for n = 1 
(Fig- IS) simply move along the VImIi axis when the assumed 
value of n is changed: 



A(f2 M /i) = -0.21A«. 



(16) 



Parameter £!/,/i 2 : This parameter affects the mass function 
model via the transfer function, but the effect is miniscule for 
reasonable values of fib /Qm- The effect of the average baryon 
density on the observed baryon mass function (via the scaling 
of Mb, eq. |4j is much stronger. This effect is equivalent to 
the /i-scaling of Mi, considered above — it shifts the erg - £ImIi 
confidence intervals along the erg axis, 



Act 8 =-18.0A1Vi 2 



(17) 



Non-universal baryon fraction, T: As we noted before, any 
large deviations of the cluster baryon fraction from the uni- 
versal value flb/H-M can significantly affect our results. 
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FIG. 8. — Constraints on <ts using only the normalization of the baryon 
mass fraction and assuming h = 0.71. The dotted lines correspond to the 
68% con fidence region and the solid line shows the analytic approximation 
(eq. 1191 1- The dashed line shows a typical <j$(£Im) dependence obtained 
from the normalization of the cluster temperature function (adapted from Sel- 
jak 2002). 



The simplest case is a uniform scaling of the baryon frac- 
tion by a factor T , i.e. T(M) = T (l -M*/M b ) for our base- 
line model (eq. 1131 or simply T(M) = To. Such a variation 
in T simply shifts the observed baryon mass fraction along 
the mass axis by a factor of Tq 5 (eq. I14> . The effect on 
the a% — QmIi confidence intervals is the same as in the case 
of variations of il^/z 2 and h — the confidence intervals shift 
along the a% axis, 

Act 8 =-0.37AT , (18) 

and the constraints on are unaffected. 

If T is a function of the cluster mass, as in dl3> . the 
shapes of the total and baryon mass functions are different 
(cf. eq.H4>. Therefore, the derived values of both erg and ft/^h 
are affected. The results of fitting the observed mass function 
assuming either T(M) given by equation ( 1131 or T(M) = 1 
are shown in Fig. [7\ For T(M) = 1, we obtain lower val- 
ues of both ilnfh and erg (dashed contour), although the effect 
on CT8 is partly due to the correlation of this parameter with 
the assumed value of £Im (see below). If the baryon fraction 
decreases in low-mass clusters, the cumulative baryon mass 
function is progressively shifted towards smaller M. This 
makes the predicted F(Mi,) flatter with a smaller normaliza- 
tion than for T = 1 . Therefore, the data will require steeper 
models for the total mass functions (i.e., higher values of 
VLmIi) with a higher normalization (i.e., higher values of <j%). 

Dotted contour in Fig. 7 shows the effect of a 10% cluster- 
to-cluster scatter in T on the og - f2 M /i constraints. As ex- 
pected, the scatter has a small effect on the results, only to 
expand the size of the confidence interval by ~ 12%. 

6.3. Constraints from Normalization of the Baryon Mass 
Function 

The <jg constraints are often derived from the cluster abun- 
dance data by using only the normalization of the mass or 
temperature functions, e.g., the number density of T > 6 keV 
clusters. We have applied this approach to our baryon mass 



Table 2 

Variations of a and b in the relation ct&(Qm) = a + bfl° M 5 as 

THE MAIN PARAMETERS ARE CHANGED. 



Parameter 


Ac, 


Ab 


h 


0.15 


0.47 


n b h 2 


-6.17 


-23.7 


T from eq. 1 1 31 


-0.08 


-0.63 


T = const 


-0.12 


-0.51 



function data. The number density of clusters with Mb > 
10 14 M Q , the median mass in the sample, was used as the only 
observable. This mass threshold approximately corresponds 
to clusters with T = 5 .5 ke V (from the - T relation pre- 

sented in Voevodkin et al. 2002). Since no information about 
the shape of the mass function is used, it is possible to con- 
strain only erg as a function of Q,m- The results are shown in 
Fig- [HI by dotted lines. The allowed band (68% confidence) 
of erg for our baseline model is accurately described by the 
following equation (solid line in Fig.[8j: 

<7 8 (f4f) = 0.60 + 0.28 f^ 5 ± 0.04. (19) 

The (78 constraints from the normalization of the cluster 
temperature or total mass functions are also il^-dependent. 
Typically, one finds crs oc fi^jj* with a = 0.4-0.6 from these 
analyses. Our f^-dependence is much weaker (solid vs. 
dashed line in Fig. [8}. The dependence of cts on in the 
baryon mass function analysis is weak, because the relation 
between the cluster baryon mass and the corresponding linear 
scale is f^-independent, A 3 ~ Mb/(pb), while for the total 
mass, A 3 ^ M tot /(flM p a )- This is one of the advantages of 
our method. The price we have to pay is the introduction of 
new dependencies for the derived value of a% on parameters 
h, flbh 2 , and T. Table |2] summarizes the changes in the coef- 
ficients of equation dl9t as these parameters are varied. 

7. COMPARISON WITH EARLIER WORK 

Determination of erg from various datasets at z = is a well- 
established field and so a comparison of our results with at 
least some recent relevant studies is in order. 

Our results can be most directly compared with the deter- 
mination of (78 from the cluster abundance data. Most such 
studies use the cluster temperature as a proxy for the total 
mass and a normalization of the local temperature function 
at T = 5 - 7 ke V as the primary observable. A review of re- 
cent results from the temperature function can be found in 
Pierpaoli et al. (2003). As we noted before, the scatter in the 
derived values of (78 in this method is mostly due to different 
assumptions regarding the normalization of the M— T relation. 
Generally, our determination of (78 is in excellent agreement 
(for f2 = 0.3) with those studies that use the M— T normaliza- 
tion from the X-ray total mass measurements assuming hydro- 
static equilibrium for the ICM (e.g., Markevitch 1998, Seljak 
2002). Our (78 is lower than the determination from the optical 
cluster mass function by Girardi et al. (1998). However, our 
value (for f2 M = 0.3) is consistent with the analyses of Bahcall 
et al. (2003) who used the cluster mass function estimated 
from the Sloan Digital Sky Survey, and of Melchiorri et al. 
(2003) who combined the SDSS results with the pre-WMAP 
CMB measurements. 

Two papers have used essentially the local X-ray luminosity 
function (XLF) and the X-ray normalized L x - M lot relation 
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(Reiprich & Btihringer 2002, Allen et al. 2003) to estimate the 
total mass function. Our results on a% are in a good agreement 
with these works (for fl = 0.3). Reiprich & Bohringer also 
fit the shape of their mass function to derive VImIi = 0.09 ± 
0.03. This is consistent but on the low side of our confidence 
intervals (Fig.[6ji, likely because Reiprich & Bohringer used 
the Press-Schechter approximation for the mass function in 
their modeling which is less accurate than the Jenkins et al. fit 
employed here. 

Matter density fluctuations can be traced by the spatial 
distribution of galaxies. Two modern surveys, SDSS and 
2dF (York et al. 2000, Colless et al. 2001) provide excellent 
datasets for these studies. The galaxy power spectrum from 
2dF implies VL M h = 0.20 ± 0.03 (Percival et al. 2001), consis- 
tent with our results for the baseline model. 

Using clusters as tracers of large scale structure is also a 
promising approach to determine the matter power spectrum. 
Using the power spectrum of 426 clusters from the REFLEX 
survey (Bohringer et al. 2001), Schuecker et al. (2003) de- 
termine a $ = 0.71 ± 0.03 and Q M h = 0.24 ± 0.02. Note that 
this method is probably less sensitive to the normalization of 
the M tot -L x relation than using the normalization of the XTF 
or XLF because the bias factor changes more slowly than the 
spatial density of clusters as a function of mass (Mo & White 
1996). It is encouraging that our results are in good agreement 
with the Schuecker et al. values. 

Finally, erg can be determined from the fluctuations of the 
cosmic microwave background. The recent WMAP measure- 
ments of the CMB angular fluctuations and Thompson optical 
depth of the Universe indicate a% = 0.9 ±0.1 assuming pure 
CDM, w = — 1, and some other common cosmological priors 
(Spergel et al. 2003); this is higher than, but consistent with, 
our value. 

8. SUMMARY AND CONCLUSIONS 

We have derived the baryon mass function for galaxy clus- 
ters at a median redshift (z) = 0.05. The baryon (gas+stars) 



mass is measured within a radius of mean baryon overdensity 
S = 324, using direct deprojection of the ROSAT X-ray imag- 
ing data (for gas mass) and an established correlation between 
M g and L opt (for stellar mass). 

We argue that if the baryon fraction in cluster is nearly uni- 
versal, the baryon mass function is an excellent proxy for 
the total mass function — F(M tot ) is obtained from Fi,(Mb) 
simply by a constant log shift along the mass axis. In fact, 
Fb(Mb) contains enough information to derive the amplitude 
and slope of the density fluctuations at cluster scales, even if 
the absolute value of M tot /M/, is unknown. 

Using this method we obtain erg = 0.72 ± 0.04 and Q,mIi = 
0.13 ±0.07. The baryon mass function method has an im- 
portant advantage in that it does not rely on rather uncertain 
observational determinations of the total cluster mass. The 
tradeoff is the sensitivity of the derived values of a% and £ImIi 
to other cosmological parameters (such as h, fi;,). Our method 
is also sensitive to any significant deviations of the cluster 
baryon fraction from universality (§|6j. However, the good 
agreement of our results with a large number of measurements 
by other, independent methods (§0 shows that our assump- 
tions about the weak trends in the baryon fraction are reason- 
able. We hope that future progress in observations and phys- 
ical models of the ICM will reduce this uncertainty to such a 
level that easily measurable baryon masses can be used as a 
reliable proxy for the total mass at all redshifts. 

We thank A. Jenkins for providing the software for com- 
puting the model mass functions, and M. Markevitch and 
C. Jones for careful reading of the manuscript. This work was 
supported by Russian Basic Reseach Foundation (grants 02- 
02-06658 and 00-02-17124) and by NASA grant NAG59217 
and contract NAS8-39073. A. Voevodkin thanks SAO for hos- 
pitality during the course of this research and also aknowl- 
edges support from the Russian Academy of Sciences via the 
Young Scientist program. 



REFERENCES 



Akritas, M., & Bershady, M. 1996, ApJ, 470, 706 
Allen, S. & Fabian, A. 1998, MNRAS, 297, L57 

Allen, S. W., Schmidt, R. W., & Fabian, A. C. 2002, MNRAS, 334, LI 1 
Arnaud, M., Neumann, D., Aghanim, N., Gastaud, R., Majerowicz, S., & 

Hurhes, J. 2001, A&A, 365, 80 
Arnaud, M., Rothenflug, R., Boulade, O., Vigroux, L. & Vangion-Flam, E. 

1992,254,49 
Bahcall, N, 1975, ApJ, 198,249 
Bahcall, N. A. et al. 2003, ApJ, 585, 182 

Beers, T., Geller, M., Huchra, J., Latham, D. & Davis, R. 1984, ApJ, 283, 33 

Bialek, J. J., Evrard, A. E., & Mohr, J. J. 2001, ApJ, 555, 597 

Blanton, E., Sarazin, C, McNamara, B., & Wise, M. 2001, 558, 15 

Bohringer, H. et al. 2001, A&A, 369, 826 

Bond. J. R. & Efstathiou, G. 1984, ApJ, 285, L45 

Buries, S., Nollett, K. M., & Turner, M. S. 2001, ApJ, 552, LI 

Colless, M. et al. 2001, MNRAS, 328, 1039 

David, L., Slyz, A., Jones, C, Forman, W., Vrtilek., S. D. & Arnaud, K. A. 

1993, ApJ, 412, 479 
Eisenstein, D. J., & Hu, W. 1998, ApJ, 496, 605 
Evrard, A. E., Metzler, C. A., & Navarro, J. F. 1996, ApJ, 469, 494 
Evrard, A. E. et al. 2002, ApJ. 573, 7 

Fabian, A., Hu, E., Cowie, L. & Grindlay, J. 1981, ApJ, 248, 47 

Feldmeier, J. J., Mihos, J. C, Morrison, H. L., Rodney, S. A., & Harding, R 

2002, ApJ, 575, 779 
Finoguenov, A., Reiprich, T. & Bohringer, H. 2001, A&A, 368, 749 
Fischer, P. & Tyson, J. 1997, ApJ, 114, 14 
Freedman, W. L. et al. 2001, ApJ, 553, 47 
Fukazawa, Y. et al. 1998, PASJ, 50, 187 

Fukugita, M., Hogan, C. J. & Peebles, P. J. E. 1998, ApJ, 503, 518 
Girardi, M., Borgani, S., Giuricin, G, Mardirossian, E, & Mezzetti, M. 1998, 
ApJ, 506, 45 

Guth, A. H. & Pi, S.-Y. 1982, Physical Review Letters, 49, 1110 
Hawking,, S. W. 1982, Phys. Letters, 115B, 295 
Henry, J. P. 2000, ApJ, 534, 565 

Horner, D. J., Mushotzky, R. F. & Scharf, C. S. 1999. ApJ, 520, 78 



Hradecky, V., Jones, C, Donnelly, R., Djorgovski, S., Gal, R. & Odewahn, S. 

2000, ApJ, 553, 521 
Hu, W. & Kravtsov, A. V., 2003, ApJ, 584, 702 

Kaastra, J., Ferrigo, C, Tamura, T., Paerels, F., Peterson, J., Mittaz, J. 2001, 
A&A, 365,99 

Jenkins, A., Frenk, C. S., White, S. D. M., Colberg, J. M., Cole, S., Evrard. 

A. E., Couchman, H. M. P., & Yoshida, N. 2001, MNRAS, 321, 372 
Lin, Y, Mohr, J. J., & Stanford, S. A. 2003, ApJ, 591, 749 
Markevitch, M. 1998, ApJ, 504, 27 
Markevitch, M. & Vikhlinin, A. 1997, ApJ, 491, 467 

Markevitch, M., Forman, W., Jones, C, Sarazin, C, & Vikhlinin, A. 1998, 
ApJ, 503, 77 

Markevitch, M., Vikhlinin, A., Forman, W., Sarazin, C, & Vikhlinin, A. 

1999, ApJ, 527, 545 
Mathiesen, B., Evrard, A. E. & Mohr, J.J 1999, ApJ, 520, L21 
Melchiorri, A., Bode, P., Bahcall, N. A., & Silk, J. 2003, ApJ, 586, LI 
Mo, H. J. & White, S. D. M. 1996, MNRAS. 282, 347 
Mohr, J. J., Mathiesen, B.. & Evrard, A. E. 1999, ApJ, 517, 627 
Mewe, R., Gronenschild, E. H. B. M. & van den Oord, G. H. J. 1985, A&AS, 

62, 197 

Nevalainen, J., Markevich, M. & Forman, W. 2000, ApJ, 532, 694 
Percival, W. J. et al. 2001, MNRAS, 327, 1297 

Peebles, P. J. E. 1980, The Large-Scale Structure of the Universe (Princeton 
Univ. Press) 

Pierpaoli, E., Borgani, S., Scott, D., & White, M. 2003, MNRAS, 342, 163 

Ponman, T., Cannon, D. & Navarro, J. 1999, Nature, 397, 135 

Pratt, G. W. & Arnaud. M. 2002, A&A, 394, 375 

Pratt, G. W. & Arnaud, M. 2003, A&A, 408, 1 

Reiprich, T. N, Bohringer, H. 2002, ApJ, 567, 716 

Seljak, U. 2002, MNRAS, 337, 769 

Seljak, U. & Zaldarriaga, M. 1996, ApJ, 469, 437 

Snowden, C, McCammon, D., Burrows, D. & Mendenhall. J. 1994, ApJ, 
424,714 

Spergel, D. N. et al. 2003, ApJS, 148, 175 
Starobinskii, A. A. 1982, Phys. Lett., 117B, 175 



BARYON MASS FUNCTION 



11 



Sun, M., et al. 2003, ApJ, in press i astro-ph/0308037 i 

Vikhlinin, A., Forman, W. & Jones, u. lyyy, ApJ, 3/5, 47 

Vikhlinin, A., VanSpeybroeck, L., Markevitch, M., Forman, W. R., & Grego, 

L. 2002, ApJ, 578, L107 
Vikhlinin, A., et al. 2003, ApJ, 590, 15 



Voevodkin, A., Vikhlinin, A. & Pavlinsky, M. 2002, Astron. Letters, 28, 793 
White, S. D. M., Navarro, J. F, Evrard, A. E., & Frenk, C. S. 1993, Nature, 
366, 429 

White, D. 2000, MNRAS, 312, 663 
York, D. G. et al. 2000, AJ, 120, 1579 



